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We study the effects of stochastic thermal fluctuations on the instability of the free surface of a 
flat liquid film upon a solid substrate. These fluctuations are represented as a standard Brownian 
motion that can be added to the deterministic equation for the film thickness within the lubrication 
approximation. Here, we consider that while the noise term is white in time, it is coloured in space. 
This allows for the introduction of a finite correlation length in the description of the randomized 
intermolecular interaction. Together with the expected spatial periodicity of the flow, we find a 
dimensionless parameter, j3, that accounts for the relative importance of the spatial correlation. We 
perform here the linear stability analysis (LSA) of the film under the Influence of both terms, and 
find the corresponding power spectra for the amplitudes of the normal modes of the instability. We 
compare this theoretical result with the numerical simulations of the complete non-linear problem, 
and find a good agreement for early times. For late times, we find that the stochastic LSA predictions 
on the dominant wavelength remains basically valid. We also use the theoretical spectra to fit 
experimental data from a nanometric melted copper film, and find the corresponding times of the 
evolution as well as the values of the parameter /3. 

PACS numbers: 


I. INTRODUCTION 


A basic problem in the study of free surfaces instabilities is the breakup of a flat thin liquid film on a solid substrate. 
Up to now, the description based on the hydrodynamic and deterministic Navier-Stokes equations are proven to be 
valid even down to the nanometric scale [T] . This has been accomplished by introducing the intermolecular interaction 
between the liquid and the substrate. However, it is known that at these scales the thermal agitation of molecules is 
relevant when describing the behavior of matter US]. Thus, it is still necessary to investigate what role can play the 
thermal fluctuations in the hydrodynamic description of these instabilities. The consequences derived from considering 
this other effect can be of interest when designing microfluidic devices or electronic components whose function relies 
on thin film properties. In particular, we are interested on the effects that thermal noise may cause on films laterally 
much larger (up to microns) than their thicknesses. 

The study of the effects of thermal noise in the hydrodynamical equations was first introduced phenomenologically 
many years ago by Landau |S] and Uhlenbeck |3]. This inclusion can be done from the deterministic Boltzmann 
equation by a long-wave approximation, which justifies its microscopic feature. These equations have been used in 
the study of turbulence in randomly stirred fluids [7], the onset of instabilities in Rayleigh-Benard convection |S] 
and Taylor-Gouette flow |3]. This subject is of interest nowadays because one of the issues to be considered in the 
discussions about the differences between Navier-Stokes equation and molecular dynamics simulations is the effect of 
thermally triggered fluctuations in the classical hydrodynamic continuum modeling. 

In particular, the application to unstable polymeric thin films has been the object of several theoretical and experi¬ 
mental studies [Mil]- In this problem, stochasticity has been analyzed using several techniques (such as Minkowsky 
invariants ini) to contrast some theoretical predictions with experimental results where stochasticity is mostly 
considered as a spatial white noise. On the other hand, the problem of unstable liquid metal films with thermal noise 
has not been object of such a thorough study. In this problem, the solid coating is melted by laser, and this introduces 
aspects that require the consideration of new factors such as the spatial correlation. Since the deposition of energy 
is not strictly uniform throughout the illuminated spot, the liquid lifetimes of different regions are not the same. In 
this context, there is a mix of factors to be considered when looking at samples from different regions. In fact, the 
time evolution of the sample at a certain region, due to the corresponding liquid lifetime, is compounded with the 
possibility that the laser illumination induces thermal fluctuations that might not be the same for all regions. One 
of the aims of the present paper is to consider how different spatial correlations of the ensuing fluctuations could 
influence the final spectra of the unstable modes. 

In this work we study the thin film instability by using a stochastic version of the thin-film equation based on 
the lubrication approximation for incompressible hydrodynamic equations (see Section 0. Then, in Section |III| 


we 


perform the linear stability analysis of the thin film under the perturbation with normal modes, and in Section IV 










2 


we solve numerically the stochastic thin-film equation, and compare the results with the linear solution obtained 
previously. A comparison of theoretical predictions with experimental Fourier spectra obtained from SEM images of 
the instability of a melted copper film is presented in Section [Vj and finally we summarize and discuss the results in 
Section |Vll 


II. THIN FILM EQUATIONS WITH STOCHASTIC NOISE 

In the framework of the continuous mechanics, the thermal agitation of the film molecules modifies the surface forces 
which describe the interaction between the fluid inside a volume element and its surroundings. Thus, an additional 
term, S, in the expression of the Newtonian stress tensor has to be considered in order to include somehow the effect 
of molecular thermal motion El US- Within the lubrication approximation, the most relevant component of S is Siz , 
where i can be either x or y and indicates a direction parallel to the substrate while z stands for the normal one. Due 
to its randomness, Siz has zero mean 


{S^z{x,t)) = 0 , 


( 1 ) 


and the correlator is given by 


{Siz{x,t)Sjz{x,t')) = 2ykBTF{x - x) 5{t - t') (2) 

where i,j = x,y, y is the fluid viscosity, and x = {x,y). Here, kg and T are the Boltzmann constant and fluid 
temperature, respectively. This property is a consequence of the fluctuation-dissipation theorem of the statistical 
mechanics, which relates the fluctuations of physical quantifies to the dissipative properties of the system. From a 
physical point of view, the hydrodynamical equations are only valid at a scale large as compared to the molecular 
scale. Strictly speaking, since thermodynamic equilibrium is characterized by a Gaussian velocity spatial distribution, 
uncorrelated noise is required. Therefore, the correlation function F{x — x') in Eq. ([^ must have a small width. In the 
same approximation, the pressure terms in the isotropic part of the stress for a film of local thickness h{x, t) are given, 
as usually, by the capillary pressure, —(where 7 is the surface tension), and the disjoining-conjoining pressure 
(van der Waals force), n(/i). Thus, the reduction of the Navier-Stokes equations under the lubrication approximation 
leads to m- 


dh 




- V 


{h - z)Snz{z)dz 


= 0 , 


( 3 ) 


where iSip = {Sxz,Syz). Note that the new noise term in Eq. ([3| enters as a rather complicated integral. It sums up 
z-uncorrelated noise terms over the film thickness, but it has fne advantage that it maintains the conservative form 
of the equation. Thus, we have now a random current which acts as another driving force. It can be shown |15j that 
the Fokker-Planck equation from Eq. (§ leads to the same time evolution of the thickness distribution function as 
that of the Langevin equation 


3/r 


dh 

m 


(yV^/i-k n(h)) 


- V 


V3fF^{x, t) 


= 0 , 


( 4 ) 


with a single multiplicative conserved noise vector ^{x,t), where the noise amplitude, ^{x,t), satishes [HI ITb] 

(|(f, f))=0, {^i{x,t) ^j{x',t')) = 2fj.kBTF{x - x') S{t - t') Sij. (5) 


In general, ^ is a correlated noise in space, but a white noise in time. 

Assuming symmetry along y-axis, the one dimensional version of Eq. Q for h{x,t) is 


„ dh d 


h^ 


f d^h 

V ^ J 


A 

dx 


V3h^^{x, t) 


= 0 , 


( 6 ) 


where, for brevity, ^{x,t) stands for ^^{x,t). 

Since the only characteristic length scale of an infinite him is its thickness, ho, we dehne the following non dimen¬ 
sional variables, 


X 



y = 7 -, ^ 

ho 


^0 ’ 



n = —n, 

7 


( 7 ) 
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where the scales of time, toi and noise, are to be determined in terms of the characteristic parameters of the problem. 
Note that we take the capillary pressure, 7 /ho, as the scale for the disjoining pressure. Thus, the dimensionless version 
of Eq. Q is: 


dh d 
di dx 


where 


d^h an 

dx'^ dx 

3fj,ho 


A r 

dx 


t) 


= 0 , 


to — 


7 


Co= ihy^. 


( 8 ) 


(9) 


The temperature scale, Tq, can now be obtained from Eq. ([^. In fact, by defining the dimensionless temperature 
T = TITq, we obtain 

(C(a;, t) i{x\ t')) = f F{x - x') S(t - t'), (10) 

with 


To = 


2kBfJ, 


( 11 ) 


Moreover, it is still convenient to define the dimensionless noise amplitude as 0 = since the correlation 

{Q{x, t) Q{x', t')) = F{x — x ) 5{t — t') (12) 

is now normalized to one. Einally, the governing Eq. ® becomes 


dh d 
di dx 


d^h an 

dx^ dx 




= 0 , 


where we define 


_ fcfiT 

2 ihr 


(13) 


(14) 


by using Eqs. and As a result, we obtain a meaningful interpretation of the dimensionless constant a. In fact, 
it gives the relative importance of the magnitude of the stochastic term (thermal noise) respect to the deterministic 
part of the equation is given by a, in the form of the ratio between the thermal and surface energies of the system. 
Since typical experimental data yield a of the order of 10“^ (or even less) we will consider here this parameter within 
this range of values in order to look for effects on the film instability. 

As regards to the form of 11 , we take into account both the attractive and repulsive intermolecular liquid-solid 
forces, so that it includes both the disjoining and conjoining pressure terms in the form 


n(/i) = Kf{h) = K 


where /i* is the dimensional equilibrium thickness, and k (with units of pressure) is given by 

A 

^ 67r/i3 


(15) 


(16) 


being A the Hamaker constant. In dimensionless variables, k becomes AT = K/ 10 / 7 , and then the final version of 
Eq. (13) is: 

dh 


dt 


d_ 

dx 




— h^^'^Q{x,t) 


dx - 


= 0 , 


where we omit the tilde (~) for brevity and for now on. 

Eor the stochastic term, we consider as usual that Q{x,t) is related to a standard Brownian motion as 

dW{xA) 


0(x,t) = 


dt 


(17) 


(18) 


which satisfies 


W{xA + Ad,-W{x,t) ~ A7(0,A), (19) 

where A/'(0, A) is a normal distribution with zero mean and variance A. Here, the notation means an equality 
of distributions. 
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III. LINEAR STABILITY ANALYSIS (LSA) OF THE STOCHASTIC THIN FILM EQUATION 


At the beginning of the instability process the deviations, 5h{x, t) = h{x, t) — hg, from the initial average film height 
are small (even if ho = 1, we keep this notation for clarity). By expanding Eq. (j^ up to first order in 6 h and 0 
(assuming that the noise amplitude is small as well) we obtain the linear stochastic equation, 


ddh -3 / d'^Sh 


KfCho) 


d'^Sh 


+ ^0 1 V-u/ I 




dt ' V 

It is convenient to look for its solution in the Fourier space, so that we have 

dSh{q, t) 
di 


I? = 0. 


= uj{q)Sh{q, t) + iy 2ahQ qQ 


where the Fourier transform is defined by 

/ OO 

6 h{x,t) dx, 

-00 

and w(g) is the deterministic dispersion relation. This one is given by |16j 

uj{q) = AuJr. 

where 


2 / \ 4 

Qc 


9c = yKf'{ho), 




( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


are the critical (marginal) wavenumber and the maximum growth rate, respectively. The wavenumber of maximum 
growth rate is qm = qc/ 

Since Eq. is an equation of the Langevin type, its solution is given by [HI US], 

Sh{q,t) = e‘^^‘^'>^Sh{q,0)+iJ2ahlq f s) . (25) 

Jo 

In order to study the instability evolution in the spectral space, we calculate the autocorrelation 

( 6 h[q, t) Sh{q', t ')) = Ai + A 2 + A 3 + A 4 , (26) 


where the terms on the r.h.s. are defined as follows: 

Ai = (^(q,0)^(g',0))e‘^(«)*e“(«')*', (27) 

A 2 oc (Sh(q,0) dW(q',t')), 

A 3 oc (Sh(q',0)dW(q,t)}, 

A4 = - 2 ahlq^{ f -^JdW{q', s')) . 

Jo Jo 

In order to calculate these terms, let us first consider the autocorrelation of the Fourier transformed noise, 0, as, 


where 


^00 /*oo 


(0(9,t)0(9',t')) = 


( 0 (a;, t) e{x', t')) ^ dx dx' 

/ OO poo 

/ 6 {t - t') F{x - x') dx dx' 

-OO j —OO 


> —OO —OO 

^00 poo 


= 27^ 5{q +q')5{t-t')F{q) 



F{u)e-"i''du 


(28) 


(29) 
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is the Fourier transform of the correlation function, F{u), being u = x — x\ and we have used e ^‘‘^dx = 27r S{q). 

Note that only in the case of non-correlated noise we have F{q) = 1, otherwise this transform has to be calculated 
(see Section III A). Since 0 is a white noise in time, its Fourier transform satisfies (see Eq. (18l) 


e{q,t) = 


dW{q,t) 
dt ■ 


Then, the autocorrelation of W{q,t) is given by 

{W{q,t)W{q\t')) = 2TT6{q + q')F{q){tAt'), 

where t A f stands for the minimum of t and t'. 

On the other hand, the height-height correlation for the initial condition is 

/ OO nOO 

/ ( Sh{x, 0) Sh{x', 0)) dx' 

-OO j —OO 

/ OO ^oo 

-OO J —OO 

= 2tt Fo{q) 5{q + q) . 

where Fq(u) is the spatial correlation function of the initial condition, and Fb(g) its Fourier transform. 
Consequently, we can write Ai in Eq. (27) as 


(30) 




where we have considered the parity uj{—q) = ui{q). 
Regarding A 2 and A 3 in Eq. (27), we note that 


Ao — A 3 — 0 


(31) 


(32) 


(33) 


because the randomness of the initial condition is independent of the Brownian W. 

Eor the term A 4 in Eq. (27) we note that since the Brownians in different time intervals t and t' are not correlated, 
only the common interval |0, t A t'] contributes to the correlation of the product of the integrals. Besides, due to 
Eq. (30) only the terms with q' = —q have non zero correlation. Thus, we obtain 


A 4 = —2a h\q^ 2 t: 5{q + q)E 


L^o 

t/\t' 




'>dW{s) 


rtAt' 


s‘^(«)(‘'-"')dIF(s') 


= -2cr 


rt/\Z 

hi q^ S{q + q')F{q) / —). 

Jo 


(34) 


The last line above is consequence of one of the lemmas of Ito’s integral [13 [E], since the stochastic process in time 
is a white noise. Performing the integral and using t + t' — 2{t A t') = \t — t'\, we have 


A4 = ahl2TTS{q + q') ^ 

uj(q) 




(35) 


Einally, by replacing Eqs. (32) 


and (35) into Eq. (^26F we obtain 


{Sh{q,t)Sh'{q',t')) = 27r (5(q + g')5'(g; C i')> 


(36) 


where 


S{q-,t,t') = Fo{q)e^^‘‘ht+t’) + 


q^F{q) 


M) 


^u>{q)(t+t') _ ^u>{q)\t-t'\ 


(37) 


Eor the case of non-correlated noise, we have F{q) = 1, which reduces this equation to that obtain in |2]- Note that 
the first term of Eq. (37) corresponds to the spectra predicted by the deterministic model (cr = 0). Here, we aim to 
compare the evolution of films with (cr > 0) and without (cr = 0) the stochastic term. To study the corresponding 
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spectra separately the film has to be perturbed at t = 0, otherwise the film does not evolve in the deterministic case. 
Thus, we consider that the originally flat free surface of the film is slightly modified by a small amplitude perturbation 
of the form, 


N 

Sh{x,0) = '^^BkSm{2Trxk/L), (38) 

fc=i 

and the amplitudes are random numbers with \Bk\ < Bmax = 10 ~^hQ. 

As a typical case, in the following calculations we choose a film with h* =0.1 and 9 = 30°, which yields [TB] qm = 
0.151, Qc = 0.213 and ijjm = 5.19 10“"^. The quantities jqm = 41.6 and = (l/wm) ln[(ho ~ /B^ax] = 

13113.5 give a rough idea of the spatial extension and time duration of the film breakup process. We find that 
L = 500 « 12Am is large enough to produce results which are practically domain size independent. The consequences 
on the stochastic process of using a correlated noise on a finite domain is analyzed in the next section. 


A. Correlated stochastic noise in a finite domain 


Here, we will assume that the correlation function F(x — x') in Eq. © is L-periodic. Moreover, we give meaning 
to the stochastic process Q{x,t) in terms of a Q-Wiener process in the form (see Eq. @) 


+ 00 


e{x,t) = ^ XkPk{t)gk{x), 

k— — oci 


dt 


(39) 


where /3fc (with k integer) forms a family of mutually independent Brownian motions with respect to time and the 
dot stands for time derivative. The constants Xk are the eigenvalues of the Hilbert-Schmidt operator Q, defined by: 


fL/2 

Qf{x) = / F{x - x')f{x')dx'. 

J-L/2 


The corresponding complete system of orthonormal eigenfunctions, gk{x), that satisfy 

Q9k{.x) = Xk9k{x) 


(40) 


(41) 


are 


gk{x) = 


1_ 

L’ 


COS (qkx) , k > 0 

k = 0 


(42) 


j sin (qkx ), fc < 0 


In fact, this can be easily verified that by considering the complex eigenfunction Gk{x) = e with qk = 27rfc/L. 
By defining u = x — x\ we obtain 


/ i/2 I'Ll j2 — X 

F{x - x')e-^^>^^'dx' = / E(u)e-*«''“dM = XkGk{x) 

-Lj2 J-L/2-X 


Lj2—x 
— Ll2—x 


(43) 


with the eigenvalue given by 


fL/2 

Xk= / F(u)e-*«'=“dM. 

J-L/2 


(44) 


Here, we used the parity property F{u) = F{—u), and the x-dependence at the limits of integration has been omitted 
due to the assumed periodicity over a distance L (which let us take x = 0 in both limits without loss of generalization). 
Thus, Eq. (441 allows to obtain all the eigenvalues for a given correlation function, F{x — x'). Note that this equation 
is the finite size domain version of Eq. ( [2^ for a discrete spectrum, so that the correlated noise effect is embedded in 
the discrete spectrum of the Hilbert-Schmidt operator Q. 




Now, we choose the particular correlation function m 
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F{u,£c) = 


f Z ^exp 


[ <5(u). 



4>0 
4 = 0 , 


(45) 


where £c is the correlation length, and Z is such that F{u,£c)du = 1. As shown in Appendix [a| we find that the 
eigenvalue in Eq.(44) becomes 


where 


F{qk) =Xk = 


L 

Wr 


Ik{a) 


(46) 


(47) 


We show in Fig. [^this eigenvalues spectrum for several values of /3 = L/(24)- Note that for /3 —>■ oo (i.e. 4 —>■ 0), we 
have Xfc —>• 1 for all k, which leads to the limiting case of white (uncorrelated) noise. For decreasing /3 (larger £cs) the 
width of the spectrum curve diminishes monotonously. The effect of the correlation region (i.e. not negligible values 
of Xk) on the film instability can be put in evidence by comparing it with the dispersion relation uj{q) as given by the 
deterministic LSA, Fq. (23) (see dashed line in Fig.[^. For /? < 4.16, all modes (stable and unstable ones) are affected 
by the noise with increasing effect on stable ones as (3 increases. On the other hand, for /? > 4.16 only unstable modes 
are affected by the thermal noise. Note that this limiting value is related to the value of £c, so that the periodicity of 
the problem, L, and the wavelength of maximum growth. Am, plays a role to determine these regions. 



FIG. 1: Linear spectrum of eigenvalues for several values of (3 obtained from Eqs. (46l-(47l. The vertical lines indicate the 
values of Qc and qm, while the dashed curve corresponds to the deterministic dispersion relation, u}{q), given by Eq. (23l. 


The actual effect of £c on the evolution of the instability is cle arly observed in the power spectrum of the perturbation, 
S{q,t), as predicted by the linear stability analysis in Section [ni| Figurej^shows S versus q at t = 200 and t = 2000 
as given by Fqs. (37) (t = t') and (38). As expected from the analysis of Fig. the inclusion of stochastic noise 
increases the amplitude of the modes q > qc (dotted vertical line) which are otherwise stable in the deterministic case. 
Note that /? = 0 (i.e. £c = oo) is coincident with this case in which cr = 0. This increment increases with /3, that is 
as the type of noise becomes closer to white noise {£c 0). Taking very large values of 4, i-e. of /3 (e.g. /? = 50), is 

equivalent to disregarding the noise at all (cr = 0), at least for short wavelengths, since both spectra are practically 
coincident for early times, and only differ at later times for smaller g’s. 

In Fig. 1^ we show the time evolution of the wavenumber of the maximum of the spectra, qmax{t), for different 
values of /3. Note that for small /3 (say (3 < 5), we find qmax ~ qm in agreement with the deterministic prediction. 
As /3 increases up to /3 « 20 we find that qmax < qm and that it approaches qm from below. For (3 > 20 the initial 
behaviour of qmax becomes closer to qm- Finally, for /? > 27, qmax > qm for all time, and it approaches qm from 
above. 





















(a)4 = 200 (b)t = 2000 


FIG. 2: Power spectrum at two different times for cr = 5 x 10“®, and several values of /3 {L — 500 and £c = 16.67,10, 5, 2.5,1.25) 
as given by the linear prediction in Eq. (371 and the initial perturbation in Eq. ( |38| |. The vertical dashed and dotted lines 
correspond to the wavenumber of maximum growth rate {qm = 2 n/Xm) and marginal stability (qc = 2ti/\c 
dotted curve shows the deterministic spectra (ct = 0 ). 


respectively. The 



FIG. 3: Power spectrum for two different values of f) for cr = 5 x 10 and several times, (a) /? = 15 (£c = 16.67), and (b) 
^ = 50 (4 = 50). 


IV. NUMERICAL IMPLEMENTATION IN A FINITE DOMAIN 


In order to understand the nonlinear effects in the film instability, we perform numerical simulations of the evolution 
of the film governed by the nonlinear Eq. (|^. The calculations are carried out in a computational domain defined by 
0 < X < L, which is divided into cells of size Ax (typically, we use Ax = 0.1 = /i* which assures convergence of the 
numerical scheme m)- 

Equation (|^ is discretized in space using a central finite difference scheme. Regarding the spatial dependence of 
the noise term, we use here only the sinusoidal modes in Eq. (42), since no flow boundary conditions are imposed at 
X = 0, L. Time discretization is performed using implicit Crank-Nicolson scheme with relaxation factor equal to 1/2. 
Thus, the time evolution of the stochastic term is performed according to Stratonovich rules. We note that all the 
results presented in this paper are fully converged, as verified by grid refinement; more details about numerical issues 
can be found in m- Due to the discretization of the equations, the minimum possible value of the correlation length 
is £c = Ax (= 0.1 in our case), since the discretized equations cannot make any correlation below this limit. 

To discretize the time-Wiener-processes in the framework of Ito-calculus, we replace Pkitn) at a time step tn by the 
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0 1000 2000 3000 4000 


t 

FIG. 4: Time evolution of the wavenumber of the maximum of the spectra, kmax, for different values of 


forward difference quotient 


Pk{tn) 


A;3 _ I5k{tn+l) - Pkjtn) 

tji 


(48) 


The difference A/3 is normal distributed and the variance is given by the time increment At„. Thus, we approximate 
Eq. m by 


A/3 _ 
Atn ~ 


(49) 


where is a computed generated random number which is approximately iV(0, l)-distributed, i.e. its histogram is 
close to a Gaussian with media zero and unity standard deviation (we used the GASDEV routine from Numerical 
Recipes m)- Altogether, the space-time discrete noise term, Eq. (p^, is given by 


1 ^ ^ 

= XkAfk9k{x), 

k- 


(50) 


where Xk is given by Eq. (46), and gk{x) by Eq. (42). Thus, Eq. (50) is used to calculate the noise term in Eq. 


Each realization of the stochastic process requires a given seed for Af. Then, some of the numerical results presented 
below correspond to a single realization, and others to the average of 20 realizations (different seeds). A typical example 
of the evolution of a film with and without noise effects for a single realization (i.e. a given seed) is shown in Fig. ^ 
Note that for the same time of the evolution, the amplitudes of the corrugations are much larger for a > 0 (Fig.[|5-h) 
than for cr = 0 (Fig.[^-d). Thus, one of the effects of the noise is to decrease the duration of the breakup process. 

In order to study how the correlated noise affects the time evolution of the instability we first concentrate on the 
time it takes to appear the first rupture of the film. By first rupture time, we mean the moment when the film first 
reaches its possible smallest value, which is h*. Figure]^ shows the time evolution of the minimum of h{x, t), namely 
hmin{t)- Glearly, as increases the breakup time, t;,, increases, such that /3 = 2.5 {^c = 100) is practically coincident 
with the case without noise {a = 0), which has the largest time. For ct > 0, this time decreases for increasing cr. 

A parameter of interest for the drop formation problem after the first breakup is the evolution of the maximum 
thickness as the final static configuration is reached. In Fig. we show the average of hmax{t) over 20 realizations 
for different values of /3. We also plot hmin(t) for reference, and define the corresponding breakup times, tf,, as 
h{tb) = 1.05/i* = 0.0105. Figure]^ shows that in fact the evolution of hmaxit) is very weakly dependent on /3 (i.e. 
ic), since the curves h^ax versus t — tk are practically superimposed. This result implies that the noise does not have 
any effect on the drop formation process after the breakup of the film, that is, during the dewetting stage following 
the pinch off. 

Now we aim to study the effects of the correlation length in both linear (early) and nonlinear (late) stages of the 
instability. To do so, we calculate the Fourier spectra of the thickness profiles for different times. In Fig. we show 
the evolution of the spectra with ic = 10 (/3 = 25) for both early and late times. All spectra correspond to an average 
of 20 realizations, and no adjusting parameter has been used (the scales for S are different to those used in previous 
sections because a different normalization was employed in the Fourier transform of the numerical results). For early 
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3- - 3- - 3- - 3 



(a)t7 = 0,t = 3000 (b)t7 = 0, i = 6000 (c)o- = 0, t = 9000 (d)o- = 0, t = 13000 



(e)o- > 0, t = 3000 (f)fT > 0, i = 6000 (g)o- > 0, t = 9000 (h)fT > 0, i = 13000 


FIG. 5: Thickness profiles, h{x,t), of a film at different times: (a)-(d) without noise (cr = 0), and (e)-(h) with noise for a 
single realization (cr = 5 x 10“®, /3 = 25, i.e. £c = 10). 



FIG. 6: (a) Maximum, hmax, and minumum, hmin, thicknesses versus time for several values of averaged over 20 realizations, 
(b) Maximum thickness, hmax, versus the shifted time t — tb, where tb = 8600, 10468, and 11122 for /3 = 25, 8.33, and 2.5 
(4 = 10, 30, and 100.) 


times, the agreement between numerics and the linear stability prediction, Eq. (37), is very good if one considers that 
some initial noise is introduced in the numerics. For larger times, the peaks of both spectra approach though the 
numerics show higher and a bit wider spectra than those predicted by LSA. A similar situation is observed for smaller 
and larger values of /3 as shown in Fig. The main difference is that LSA overestimate the amplitude of the peaks 
respect to the numerical ones for small P (Fig. |^), but the contrary occurs for large /3 (Fig.|^). 


V. COMPARISON WITH EXPERIMENTS 

Previous comparisons between experiments and stochastic models have studied the instability of polymeric films 
on silicon oxide substrates [121 HI! ■ However, these comparisons were made without considering spatial correlation, 
i.e. assuming both spatial and temporal white noise. Also, they mainly employed the integration of the spectra S{q) 
for all possible values of g, and derived quantities from it. Here, instead, we apply the theoretical model described 
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(a)/3 = 25, early times (t>)/3 = 25, late times 

FIG. 7: Numerical power spectra, S{q,t) (solid lines), for (a) early and (b) late times for ct = 5 x 10“® and Ic = 10 (/3 = 25) 
averaged for 20 realizations of the problem defined in Fig. [^. The dashed lines are the corresponding predictions of the LSA. 
The vertical dashed line corresponds to the wavenumber of maximum growth in the deterministic case, qm = 0.151, while the 
dotted one to the marginal value, qc = 0.215. 



FIG. 8: Numerical power spectra, S{q,t) (solid lines), for ct = 5 x 10“® and (a)£c = 30 (/3 = 8.33) and (b) £c = 2 {13 = 125) 
averaged for 20 realizations. The dashed lines are the corresponding predictions of the LSA. The vertical dashed line corresponds 
to the wavenumber of maximum growth in the deterministic case, qm ~ 0.151, while the dotted one to the marginal value, 
qc = 0.215. 


above to experimental results for unstable liquid metal films to evaluate the importance of spatial correlations when 
considering stochastic instabilities. In order to do this, we do not restrict ourselves to some integrals of the spectra, 
but employ their complete profiles as a function of the wavenumber, q. 

Our experimental data correspond to copper thin films of a few nanometers thick that are melted by the illumination 
with pulses of an Excimer laser that last some tens of nanoseconds. During these pulses, the metal is in a liquid state, 
and thus the present hydrodynamic model can be applied. In this configuration, the liquid lifetime of the melted 
copper is related with the local temperature of the film, i.e. with the spatial distribution of the laser intensity, which 
spans in a radially symmetric Gaussian profile. After the pulse, the metal solidifies leaving a distinct pattern of holes, 
drops and/or ridges depending on how long the metal has been in the liquid state. More information about this setup 
configuration and details on the technique can be found elsewhere [22h26j . 

Since the outer regions of the laser spot have shorter liquid lifetimes, one can associate these regions with earlier 
times of the evolution, and consequently, central regions with later times. Since the laser spot is relatively large, 
the SEM images of these experiments have the advantage of offering more spatial information than other setups [H] . 
Nevertheless, they have the drawback that the times corresponding to every stage of the evolution are unknown, even 
if it is possible to order the time sequence in connection with the distance of the image respect to the center of the 
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laser spot m- The goal of the following comparison is to show that the experimental observations represented by 
the spectra require not only a stochastic temporal evolution, but also some spatial correlation in the thermal noise in 
order to reproduce the full results. 

In particular, we will concentrate here on the data reported in where the SEM images of the evolving melted 
metal were analyzed by using bidimensional (2D) discrete Fourier transform (DFT). Since, the 2D spectra turned out 
to be radially symmetric in the wavenumber space, {qx, Qy), the results in Fig. 5 of [27] were reported as amplitudes 
^ 20 : versus k = {q^ + These amplitudes were in fact averaged on circles of radius k, and therefore the 

corresponding ID amplitude is obtained as Aid = (see symbols in Fig. |^. The symbols for both small k and 

amplitudes {S < 0.15) are an artifact of the finite length of the sample in the Fourier calculation. 






FIG. 9: Experimental power spectra, A\D{k,t), (symbols) from Fig. 5 of [27|, and theoretical spectra (solid lines) obtained 
with the present stochastic model with spatial correlation. The experimental spectra are organized in decreasing order of their 
distance to the center of the laser spot. 

The parameters for liquid copper are 7 = 1.304 N/m, and n = 4.38 mPa s. Assuming T = 1500 K as a typical 
temperature of the film with thickness Hq = 8 nm, we have a = 2.48 x 10“"^, and to = 0.08 ns. Regarding the 
intermolecular interaction with Si02 we use {n,m) = (3,2), h* = 0.1 nm and A = 2.58 x 10“^®J (as suggested 
in [22]). Thus, we have qc = 63.4 and qm = 44.8 (dotted and dashed lines in Fig. 

In order to perform the comparison of the experimental and theoretical spectra (see Eq. (|^ I) we choose a constant 
value for the unknown Fo{q), namely Fo{q) = 2 x 10“"^, and use the same normalization factor for the DFT as in [27j . 
Thus, we are left only with t and /3 as adjustable parameters. The fitting values for the spectra in Fig. [^ are given 
in Table [U The low local maximum for k « 100/j,to“^ is related to the size of the drops, which is smaller than the 
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distance between them m- 

Interestingly, we find not only increasing values of time as one moves from inner to outer regions (as expected), 
but a decrease of the corresponding values of (3 is also required for the fitting. This implies that the stochastic noise 
is somehow different at the sampled regions which, in turn, correspond to distinct liquid lifetimes. However, the 
relatively large values of (3 for the first three images suggest that the noise is practically white at the beginning, and 
that spatial correlation becomes important only for larger times when j3 decreases significantly. In general, it is then 
expected that the spectrum for earlier times (i.e., near the outer borders of the laser spot) correspond to a quasi 
white noise, but the noise becomes more and more spatially correlated as one goes to the center of the spot (i.e. as 
the liquid lifetimes increase). In fact, the correlation length, £c, can be estimated considering the value of (3 and the 
length of the image, which can be assumed as the periodicity length, L. For the images corresponding to Fig. |^we 
have L = 2.13 /iw, so that we obtain = L/(2j3) as shown in Table [ij Moreover, note that finally approaches 
(= 144nm), which is also very close to (= 165nm). Thus, Ic turns out to be very close to the average distance 
between drops. 


Fig. 

t(ns) 

3 

tc (nm) 


a;;?’’ (nm 

) A“74 

(a) 

0.08 

175 

6.1 

22.9 

62.8 

10.3 

(b) 

0.48 

160 

6.6 

21.0 

99.7 

14.9 

(c) 

0.97 

140 

7.6 

18.4 

125.6 

16.5 

(d) 

6.21 

7.4 

144.1 

0.97 

165.3 

1.15 


TABLE I: Best fit values from the comparison of the stochastic model with spatial correlation with experimental spectra of 
unstable liquid metal films. Here, we have Xm = 144 nm. 


VI. SUMMARY AND CONCLUSIONS 

In this work we have considered the effect of correlated thermal noise on the instability of a liquid thin film under 
the action of viscous, capillary and intermolecular forces by adding a stochastic term in the lubrication approximation 
equation for the film thickness. This term depends on the noise amplitude that is spatially self-correlated within a 
characteristic microscopic distance, ic- The linear stability analysis (LSA) of the resulting equation shows that this 
yields a new factor in the stochastic part of the instability spectrum, which is given by the Fourier transform of the 
correlation function that can be expressed in terms of the eigenvalues of the Hilbert operator associated with it. 

In order to observe the nonlinear effects on the evolution of the instability, we also perform numerical simulations 
of the full lubrication equation using different seeds to generate the random sequence of amplitudes for the stochastic 
term (so that a realization corresponds to each seed), and average the resulting power spectra to obtain a representative 
spectrum to be compared with the one predicted by the LSA. As expected, we find a good agreement with LSA for early 
times. Interestingly, for late times we obtain that the wavenumber of the maximum of the spectra tends to approach 
the deterministic value, Qm, corresponding to the LSA without stochasticity. Since the LSA with stochasticity also 
tends to we can conclude that the typical lengths of the patterns in advanced stages of the instability with 
stochasticity seem to be close to the length of maximum growth rate of the linear deterministic modes. 

Therefore, encouraged by this result we also compare the LSA prediction with the experimental data from the 
instability of melted copper films on a silicon oxide substrate. These data correspond to the early stages, where the 
holes start to grow, as well as to the stages of drops formation, i.e. after having passed through the processes of 
film breakup and dewetting. A special feature of these data is that they come from different spatial regions of the 
laser spot, and thus received distinct illuminations. Thus, different times of a single evolution can be attributed to 
each region. These times were estimated here by fitting the LSA power spectrum to each experimental one with its 
corresponding value of j3. As a result, we found that the early stages of this experiment evolved with a practically 
white noise in space, while a strong spatial correlation appeared in the spectra for late times. This shows that 
the explanation of experimental results in the nanometric scale requires the inclusion of some thermal noise in the 
modeling. In particular, correlated noise seems to be an important factor in the central regions of the laser spot, i.e. 
those with larger liquid lifetimes. We believe that our results justify further testing with more detailed experimental 
data. 
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Appendix A: Eigenvalues of the correlation function 


Here, we calculate the eigenvalues of the Hilbert-Schmidt operator Q as defined by given by Eqs. (44) and (45). By 
using the variable v = 'ku/L, the eigenvalues can be written as 


Xk = 


A(a, k) 

A(a,0) 


where 


A{a,k)= [ 

Jo 


(Al) 


(A2) 


and a is given by Eq. (47). In order to perform the above integral, we make the change of variables 2v = 9 which 
leads to the following expression, 

sin^ V = -(1 — cos2n) = -(1 + sin0). 


This one allows us to write Eq. (A2) in terms of sin0, as 


A(a, fc) = ^ 


H)^-c 




The above substitution is convenient in view of the the relation, 

OO 

e“sine ^ 

— OO 

which becomes useful here upon defining x = ia. Thus, we have 

OO 

— OO 

where is the Bessel function of order m. Now, we can also use the property 


(A3) 


(A4) 


(A5) 


(A6) 


where Imia) is the modified Bessel function of order m. By replacing Eqs. (A5) and (A6) into Eq. (A3), we obtain 

/_ „\k °° 

A(a,fc) = (A7) 

— OO ~ 2 

Since, the above integral yields 2Tr5km, we finally have 

A{a,k) = TTe~°‘Ik{a), (A8) 


so that the eigenvalue in Eq.(Al) becomes 


xiQk) =Xk = 


4(a) 
4 (a)’ 


(A9) 


which is the expression in Eq. (46). 
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